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Abstract 

The disorder driven quantum Hall to insulator transition is investigated for a two-dimensional lattice system. 
We consider a Gaussian correlated random potential, study the behaviour of the current carrying states and 
trace their energetical position when the disorder strength is increased. Our results qualitatively resemble 
those obtained previously for exponentially correlated disorder potentials. We find both the downward 
movement of the anti-Chern states as well as the floating up of the Chern states across the Landau gap 
which is sometimes masked by the global broadening of the tight binding band. 
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1. Introduction 

The disorder driven quantum Hall to insu- 
lator transition once again became a topic of 
active research recently. The theoretical activ- 
ities undertaken over the past years resulted in 
a controversial debate concerning the fate of 
the current carrying states [|1], |, g, §, 0]. 
Likewise, the experimental results reported all 
along remain inconclusive [T^ |TT], [T^ 
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T^ ]. In numerical studies using un- 
corrected random disorder potentials, it was 
found that in contrast to the continuum model, 
where the critical states are known to float 
up in energy [ p^ , |17|, [1^], the current carry- 
ing states are annihilated by anti-Chern states 
originating from the band centre in the lat- 
tice model [0, ^ ^. It has also been claimed 
^ that such a scenario should be able to 
explain the direct QH to insulator transitions 
from higher Hall plateaus with filling factor 
z/ > 2, which apparently have been observed 
in several experiments. However, according to 



the proposed global phase diagram [[191], which 
is based on the levitation picture, such transi- 
tions should not be allowed. 



In a recent paper we have shown that 
the lattice model can exhibit similar proper- 
ties as the continuum model if long range cor- 
related random disorder potentials are consid- 
ered. Using exponentially correlated random 
potentials to model the intrinsic scattering po- 
tentials and other imperfections that may in- 
fluence the movement of the charge carriers in 
real 2DEGs, we found that the floating up in 
energy can already be observed for correlation 
lengths larger than about half the lattice con- 
stant. In particular, with increasing disorder 
the floating up of the critical energy across the 
Landau gap without merging has been demon- 
strated. The latter is in contrast to the results 
of a numerical study [ pi] ] where Gaussian cor- 
related disorder potentials were used. In the 
present paper we investigate whether the ob- 
served discrepancy is due to the differences in 
the correlated disorder potentials considered. 
Therefore, in the following we also use Gaus- 
sian correlated disorder potentials and calcu- 
late numerically the disorder dependence of 
the density of states and investigate the finite 
size scaling of the localisation length from 
which the energetical position of the current 
carrying states can be extracted. 
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2. Model and Method 

The QHE system is described by a dis- 
ordered two-dimensional single band lattice 
model H = E 

"Yli <mn>^n'nT-^m^n ~^ ^mnCnCj^), WherC C^, 

are the creation and annihilation operators of 
a particle at site m, respectively. The perpen- 
dicular magnetic field is incorporated via the 
complex phases of the transfer terms, Vmn = 
V exp(i6mn), between neighbouring sites m, n 
on the lattice. In the Landau gauge, 6mn = 
±27r(p/g)(f„ ■ ey)/a,iirn = rm± ae^, and 
bmn = else, where fm is the position of site 
m and e^, Cy are unit vectors pointing in the 
X and y directions. We choose p/q = 1/8 
flux quanta h/e per plaquette a x a, where 
a is the lattice constant. The Gaussian cor- 
related disorder potentials Wm are generated 
from uncorrected random variables associ- 
ated with each lattice site n by local averaging 
using a Gaussian weighting function, Wm = 



1/^ En^nexp( 



n\ 



Y?7^), with corre- 



lation length f] and some normalisation factor 
A^. This leads to a spatial correlation func- 
tion decaying qualitatively like {wmWn) ~ 
exp(— |m — n\'^ /2rf). The probability density 
of the En is taken to be P{en) = 1/{2W) in 
the range —W < < W, and zero else. 

The density of states p{E, W) is calculated 
from the distribution of energy eigenvalues 
which are obtained by direct diagonalisation 
of square systems averaged over 60 realisa- 
tions of the correlated disorder potentials. The 
exponential localisation length \jl{E, W) = 
-limi^oo(2i:)-MnTr{|Gii|2} is calculated 
numerically by means of a recursive Green 
function method [g3|,0. G = {E-H+i-f)~^ 
is the total Green function and Gil is defined 
on the subspace of the first and last columns 
on the lattice. The system's length and width 
are L and M, respectively. Both, the den- 
sity of states and the localisation length were 
computed for various system sizes, disorder 
strengths, and correlation lengths of the dis- 




Fig. 1. The density of states p{E, W) versus energy for 
disorder strengths M^/F = 1.0, 1.5, 2.5, and 3.5. The 
system size is Al/a — 48 and the lines simply connect 
the calculated data points. 



order potentials. In what follows, we restrict 
the presentation of our results to those data ob- 
tained for correlation parameter r]/a = 1.0. 



The same value has been utilised in Ref. [21]. 



3. Results and Discussion 

The movement of the energetic position of 
the current carrying states induced by a change 
of the disorder strength is obscured by the 
global broadening of the tight binding band 
with increasing disorder. To distinguish this 
effect, which in turn shifts the Landau bands 
outwards, from the true floating of the criti- 
cal states, the density of states p{E) is calcu- 
lated first for several disorder strength in the 
range 0.5 < W/V < 5.0. As an example, the 
result of p{E) for disorder strength W/V = 
1.0, 1.5, 2.5, and 3. 5is plotted in Fig. [1]. Due to 
the symmetry with respect to E/V = 0, only 
the lower half of the total band is shown. With 
increasing disorder the sub-bands broaden and 
shift slightly outwards. The effect is stronger 
near the band edge. Similar curves were ob- 
tained for further disorder strengths. The po- 
sition of the sub-band peaks were extracted 
from a curve-fitting procedure using Gaussian 
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Fig. 2. The energy dependence of the normalised lo- 
caHsation length Xm{E, W)/M for disorder strengths 
W/V = 1.5, 2.5, 3.5, 4.0, 5.5, 6.0, and 6.5. The system 
length is at least i/a = 3- 10'' and the width M /a 48. 

or power-like shapes for the particular Landau 
bands, exploiting the normalisation condition 
for each single band and tracing the fit param- 
eters from lower to higher disorder. 

The energy dependence of the normalised 
localisation length \MiE,W) /M has been 
calculated for disorder strengths in the range 
1.0 < W/V < 6.5. Curves for W/V = 1.5, 
2.5, 3.5, 4.0, 5.5, 6.0, and 6.5 are plotted in 
Fig. ^ whereby the data points shown are aver- 
ages over 30 disorder realizations. The system 
width is M/ a = 48 and the length varies be- 
tween L/a = 3 ■ 10^ and L/a = 5-10^ depend- 
ing on the varying convergence for different 
energy and disorder strength. For fixed sys- 
tem width M, Xm{E,W)/M increases with 
increasing disorder strength. Associated with 
this increase is a shift of the peak position, 
which corresponds to the floating up of the 
current carrying state to higher energies. 

In Fig. 1^ the energy and disorder depen- 
dence of p{E, W) and Xm{E, W)/M are com- 
bined. The position of the sub-band peaks of 
the density of states (x) is shown within the 
disorder-energy plane. With increasing disor- 
der, the density of state peaks shift to lower 
energy. This movement is well fitted by the 
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Fig. 3. The position of the sub-band peaks in the den- 
sity of states p{E, W) (x) and the position of the nor- 
malised locaUsation length \m{E^W)IM (O) within 
the disorder-energy plain. The anti-Chern states are 
depicted by (•). The dotted line is an empirical fit 
Wp{E) = A{En - Ey/'^ to the peak position of the 
sub-bands. 

empirical relation Wp{E) = A{En - Ef/'^ 
which corresponds to the observation that the 
bandwidth of the total band increases with dis- 
order ~ W^. The En denote the energies of 
the sub-bands (Landau levels) for W/V = 0. 
In addition to the sub-band peaks, the posi- 
tion of the current carrying states (o) are plot- 
ted also in Fig. 0. The energies of the critical 
states were extracted from a finite size scal- 
ing analysis combined with a similar curve- 
fitting procedure as used for the total density 
of states above. For each disorder strength 
\m{E,W) /M decreases for localised states 
with increasing system size while it stays con- 
stant in the thermodynamic limit only at those 
energies that correspond to the current carry- 
ing states. 

For small disorder W both the energy of the 
peaks in the density of states and the energy 
of the critical states coincide. However, with 
increasing disorder strength, the curves move 
apart. The floating of the energy of the criti- 
cal states in the opposite direction as the sub- 
bands, almost across the lowest Landau gaps. 
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is clearly seen. The downward movement of 
the anti-Chern state (•) can easily be followed 
until it is finally annihilated by the Chem state 
of the lowest sub-band at a disorder strength 
of W/V ~ 7. This also defines the cross-over 
from the quantum Hall liquid state to the insu- 
lator. 

Comparing our results with those published 
previously for exponentially correlated disor- 
der potentials [ p4l , we notice that the be- 
haviour of the current carrying states is qual- 
itatively similar. With increasing correlations 
the position where the anti-Chem and the last 
Chem state meet is shifted to higher energies. 
Therefore, it is to be expected that for stronger 
correlations only the floating up of the energy 
of the current carrying states is observed as 
in the continuum model. The apparent merg- 
ing and subsequent joint floating of the cur- 
rent carrying states as reported pT| ] for simi- 
lar Gaussian correlated disorder potentials was 
not observed in our investigation. However, 
this may look different for smaller magnetic 
field as applied in pT| ] because the distinction 
of clustered critical states becomes increas- 
ingly problematic in finite size studies. 

In conclusion, the disorder driven quantum 
Hall to insulator transition in a lattice model 
with Gaussian correlated random potentials 
has been investigated. For a correlation length 
7] /a = 1 we have shown that the energy of 
the lowest current carrying state floats up in 
energy across the Landau gap when the disor- 
der strength is increased. This is similar to the 
levitation scenario proposed for the continuum 
model. In view of these results, the proposed 
explanation pT| ] for the apparent direct transi- 
tions to the insulator from plateaus with u > 2 
as being due to the downward movement of the 
anti-Chern states is unlikely. We believe that 
the direct transitions reported in experiments 
are either unresolved single one by one transi- 
tions or have to be accounted for by electron- 
electron interactions. 
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